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A  GENERAL  INVERSE  SOLUTION  FOR  2-D  EDGE  WAVES 

K.P.  Shaw,  Professor,  S.U.N.Y.  at  Buffalo 
Y.K.  Sun,  Research  Assistant,  S.U.N.Y.  at  Buffalo 


Abstract 

An  inverse  method  is  used  to  determine  topographies  which  will 
support  long  period  edge  waves  along  a  straight  coastline.  Two  examples  are 
given;  one  corresponds  to  a  known  solution  and  the  other  appears  to  be  new. 


I :  Formulation  and  Solution 

Consider  the  conditions  under  which  a  long  wavelength  free  surface 
gravity  wave  may  be  trapped  against  a  straight  coastline  by  a  topography  which 
varies  in  a  direction  normal  to  the  coastline.  Taking  x  in  this  normal 
direction  and  v  along  the  coastline,  the  linearized  shallow  water-long 
wavelength  equations  for  a  homogeneous  fluid  on  a  rotating  earth  are 
[e.g.  LeBlond  and  Mysak  (197b)]: 

lu/lt  +  g  3n/0x  =  iv 

iv/Ot  +  g  jq/dy  =  -fu  (1) 

jq/3t  +  H(3u/dx  +  3v/3y)  +  u3H/3x  =  0 

where  (u,v)  are  velocities,  o  is  the  free  surface  elevation,  H  is  the  water 
depth  and  f  is  the  Coriolis  parameter. 

Assume  solutions  which  propagate  along  th-  coastline,  e.g. 
u(x,v,t)  =  U(x)  cos  (ky  -  at) 

v(x,y,t)  =  V (x)  sin  (ky  -  ot) 

n(x,y,t)  =  E(x)  sin  (ky  -  ot) 


Clearly  li  and  V  may  be  eliminated,  leaving  a  single  equation  on  E, 


(  2  ) 


e.g.  LeBlond  and  Mysak  (1978  -  p.  220) 
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This  equation  does  not  force  a  particular  x  dependence  on  either  E  or  H,  as  is 
the  case  for  the  three  dimensional  problem,  e.g.  Shaw  (1979).  If,  for  example, 
E  is  required  to  be  exponentially  decaying  such  as  E^exp (-£x) ,  equation  (2) 


requires  H(x)  to  be 


az  -  f 


H(x)  =  gTk2-T2)-  +  Hoexp 


(k2-£2)  1  1 

£  +  fk/oJXJ 


where  is  an  arbitrary  constant  defined  in  general  either  by  choosing 

H(0)  =  0  or  U(0)  =  0  to  have  no  net  mass  flux  across  the  coastline,  x  =  0. 
If  H(0)  is  chosen  to  be  zero,  the  topography  is  given  by 


H(x) 


a2  -  f2 

g(k?— £2) 


exp 


(3) 


corresponding  to  the  lowest  (n=0)  mode  given  by  Ball  (1967). 

Certain  ranges  of  parameters  give  negative  topographies  which  must  be  excluded. 
Requiring  that  £  >  0,  but  recognizing  that  k  may  be  either  positive  or  negative, 
yields  the  following  results: 

For  a  >  f,  k  >  £  gives  a  positive  concave  upwards  topography, 
exponentially  approaching  a  uniform  depth,  H  equal  to  (a2-f2)/g(k 2-£2) ; 

CO 

-£  <  k  <  £  gives  a  positive,  convex  upwards,  exponentially  increasing  depth, 
-o£/f  <  k  <  -£  gives  a  positive,  concave  upwards  topography  again  exponentially 
asymptotic  to  and  finally  k  <  -o£/f  yields  a  negative  topography. 

For  a  <  f,  k  >  £  and  -o£/f  <  k  <  £  both  yield  negative  topographies, 
but  -£  <  k  <  -o£/f  yields  a  positive,  concave  upwards  topography  again  decaying 
exponentially  to  the  asymptotic  depth  H  while  k  <  -£  yields  a  positive, 
convex  upwards  exponentially  increasing  depth. 


The  choice  of  U(0)  equal  to  zero  requires  either  o  =  -fk/i,  leading 
to  a  constant  depth,  H  =  f2/g£2,  or  the  trivial  solution,  E(0)  equal  to  zero 
and  thus  E(x)  equal  to  zero  everywhere.  This  choice  is  then  rejected  for 
this  particular  form  for  E(x). 


Identifying  the  asymptotic  depth  as 

H co  =  (a2  -  f2)/g(k2  -  l2)  (4) 

and  using  the  slope  at  the  coastline,  as 

=  dH(x)/dx|x=o  =  H^(k2  -  12)/(1  +  fk/o)  (5) 

provides  a  topography 

H(x)  -  Hjl  -  exp (-H^  x/lO)  (6) 

to  support  an  edge  wave 

E(x)  =  E0  exp  (-lx).  (7) 


The  choice  of  the  two  parameters  H,*,  and  \VQ  completely  fixes  the 
topography;  the  remaining  four  parameters  (a,  k,  1  and  f)  are  then  related  by 
two  equations,  i.e.  eqs.  (4)  and  (5),  leaving  two  parameters  free.  If  the 
special  case  of  no  rotation,  f=0,  is  considered,  one  degree  of  freedom  remains 
and  these  equations  may  be  written  as  a  dispersion  equation 

oV(gH^)2  +  o2/  gH<„  =  k2  (8) 

and  an  equation  governing  the  decay  parameter 

22  =  k2  "  a 2 / 8 Hco  (9) 

Thus,  an  arbitrary  frequency  a  determines  specific  values  for  k  and  Z  on  a 
given  topography.  These  match  equations  given  by  Ball  (1967).  The  important 
point  here  is  that  the  same  topography  exists  for  all  frequencies;  only  the 
trapped  wave  changes  as  a  changes. 
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It  is  of  interest  to  note  that  the  solution  due  to  Ball  (n=Q  case)  rep¬ 
resents  only  part  of  the  solution  obtained  here.  Taking  the  case  of  no 
rotation  for  simplicity,  the  Ball  (n=0  case)  represents  the  solution  for 
k2  >  £2  which  has  a  topography  which  decays  exponentially  to  H^.  The 

solution  for  k2  <  l2  given  here  has  H  <  0  and  exp{-(k2-£2)x/£}  >,  1  such 

00 

that  the  topography  is  exponentially  increasing  with  x  ,  Shaw  (1981). 

In  general,  equation  (2)  may  be  considered  to  define  E  for  a  given  H 
or  alternatively  H  for  a  given  E.  While  the  first  view  is  more  realistic 
physically,  i.e.  a  direct  problem,  the  second  is  simpler  mathematically,  i.e. 


(10) 


may  be  solved  completely  for  H(x)  for  •  -iven  E. 

.  r  r  E"(x)-k2E(x>  ,i 


H(x)  =  C^exp  - 

L 

.  if2-c2l 


E'  (x)-—  E(x) 

G 


'f2-c2]  r  [  E'  '  (x)j^k_2E£x) _  | 

+  — r-_expj-  l  :-k  dx 

<•  s  J  ^  J  e’  (x)-~-  E(x) 


I  - -  *  exp 

J  E*(x)  -  ~  E(x) 


"  f  E' ' (x)-k2E(x) 

T  FlT 


dx  dx 


where  Cq  is  an  arbitrary  constant  to  be  chosen  to  have  zero  net 
across  the  coastline,  which  is  not  necessarily  chosen  to  be  x^O,  or,  as  will 
be  seen  in  the  example,  to  keep  finite  depths  at  the  critical  point  of  equa¬ 
tion  (10).  While  eq.  (11)  provides  a  complete  solution  to  the  inverse  problem, 
it  is  more  instructire  to  integrate  eq.  (10)  directly  for  any  given  wave  form, 
E(x),  especially  when  dealing  with  the  critical  point  at  which  E'  -  (fk/o)E 
vanishes. 

Consider  as  an  example  the  form 


E(x)  ■  E_Jf(x+b)  exp[-  J(x+a)  ] 


(12) 


which  leads  to  an  analytical  solution  lor  H(x)  in  terms  of  either  an  incom¬ 
plete  gamma  function  or  a  Kummer  function.  A  nondimens ional  form  using 

_  y? 

E  -  E/EjexpUCb-a)],  H  “  H/Hq,x  =  2(x-rb),  f  -  f/c,  k  =  k /l  and  o  =  (o)/A( gHQ)  , 
and  then  dropping  the  bars  for  convenience,  gives  eq.  (5)  in  the  form 

(13) 

(l-(l+fk)x)H'  +  (-2+(l-k^)x)H  =  .;2(f2-l)x 


The  comp  lumen  Cary  solution  lo  this  equation  is 

-  01 

H  (x)  -  C  { j (l+fk)x-l I *expl~(k2-l)x/ (1+fk) ] f 
c 


(14) 


where  a  is  (1  +  2fk  +  k2)/(l  +  2l‘k  +  f2 k^).  Variation  of  parameters,  using 


G  as  G(x),  l^ads  to 


dG 

dx 


« 

-o'  ( 1  —  f  2  )  x  [  1—  ( 1  +  fk;>:J  *oxp  [  (k^-l)x/  (1  +  fk)2] 


(13) 


The  integration  i or  this  is  best  carried  out  separately  for  x  less  than  and 
greater  than  the  critical  value,  (1  +  rk)"1.  Requiring  that  H  remain  finite 
at  this  critical  value  fixes  the  constant  of  integration  for  G(x)  and  the  same 
result  is  obtained  in  both  regions 


d(x) 


•  (1-f2) 


2r(  0 

1  +  Ik 


jxp (a ) y* (a, 1 ) [ 


(16) 


(k-  -  1)  i 

where  X  =  (k‘:-l)(i-(l  +  fk)x)/(l  +  fk)2  and  is  cliat  form  of  the 

incomplete  gamma  function  which  is  analytic  for  all  A;  e.g.  Abramuwitz  and 
Stegun  (1964) 


7  *  (<X ,  '  )  =  X 


-A 


re*) 


exp(-t)tl  1  dt 


(17) 


which  is  especially  useful  for  negative  values  of  a.  in  particular,  y*(a,A) 

is  real  for  negative  values  of  \,i.c.  for  x  greater  than  the  critical  value. 

For  some  ca.-es,  nix)  may  be  shown  to  b<-  monotonic  and  asymptotic  to 

Ho  =  o2(]-f')/(k2-l)  as  x  becomes  infinitely  large.  The  origin  for  x  is 

arbi  trarv;  roa.»t  i  i  r.e  is  chosen  to  nave  hL\ 

i  coast  * 


i-iiual  Lo  zero- 
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Kor  computations,  a  series  valid  for  all  !z|  •'  >"  is  given  by  Abramowitz  and 
Stegun  (1984): 


Y* (a, z)  =  exp(-z)  l  z  /T(a  +  n  +  1) 
n=0 

CO 

=  (r(a))-1  y  (-z)n/(a+n)n! 
n-0 


(18-a) 


(18-b) 


Y*(a,z)  may  also  be  written  in  terms  of  a  Kummer  function  as 


Y*  (a,z)  =  M(l,l+a,z)/  P  C  <X+  exf(^) 


(18-c) 


Returning  to  dimensional  variables,  using  £  =  A  as  the  independent 

spatial  variable,  leads  to  a  topography 


H(0  =  H  {1  -  2r(a)(l+v)"1exp(4/HjY*(a,5/H  )> 


(19-a) 


H  U  -  2a_1(l+v)_1M(l,l+a,C/H  )} 


(19-b) 


where  v  =  (f/o)(k/£),  a  =  (1  +  2v  +  (k/£)z) /(1+v)2  and  =  (o2-f2) /g(k2-£2) 
with  all  values,  e.g.  f,  a,  k,  l,  dimensional.  The  trapped  wave  is 


where 


E(x)  =  E^(x+b)exp[-Jt(x+a)  ] 


x  =  f^l+v)'1  -  b  -  (l+v)«.C/Htt>(k2-i2) 


In  this  case,  the  topography  is  fixed  by  the  choice  of  three  parameters, 
H^,  v  and  a,  at  least  as  a  function  of  the  variable  £.  The  four  parameters  o, 
k,  4,  and  f  are  then  related  by  three  equations 


v  =  fk/oi. 

a  =  (1  +  2v  +  (k/4) 2) / (1+v) 2 
H  =  (a2-f2) /g(k2-42) 


leaving  one  degree  of  freedom,  e.g.  o  may  be  specified  arbitrarily  and 


k,  £  and  f  then  determined  from  these  three  equations  for  a  fixed  topography 


given  by  specific  H^,  v  and  a.  This  is  not  as  convenient  as  the  previous  solu¬ 
tion  which  had  two  degrees  of  freedom;  generally  f  is  determined  by  the  latitude 
of  the  observed  wave.  If  a  is  specified  here,  f  will  be  determined  by  the 
topography  rather  than  the  geography.  The  fact  that  the  topography  is  given 
in  terms  of  5  simply  means  that  the  form  of  the  edge  wave  must  also  be  expres¬ 
sed  in  this  variable,  i.e. 


E(5)  =  e  1  (1+v)  1- (l+vHcnJck2— Jl2)-1) 

*exp(-£(a-b+£-1 (1+v)"j-(1+v)£^h“1 (k2-£2)-1)  (25) 

It  is  of  interest  to  note  that  the  location  of  the  coastline,  H(£c)  =  0, 
gives  £  /H  as  a  function  of  a  and  v  alone; 

C  OQ 


cxp(£  /H  )y*(o,5  JVJ  =  (l+v)/2r(a) 

(26-a) 

M(l,l+a,£  /H  )  =  a (1+v) /2 
c  00 

(26-b) 

The  slope  at  the  coastline,  ,  as  given  by 


1  _ 

rH  'l 

f  2 

Uc) 

l  1+vJ 

-  1 


(27) 


is  also  a  function  of  a  and  v. 

Again,  for  the  special  case  of  no  rotation,  f=0,  the  equations  simplify; 
since  the  topography  depends  now  on  only  two  parameters,  and  a,  the  three 
parameters  o,  k  and  £  still  possess  one  degree  of  freedom,  i.e.  for  a  given 
and  a,  and  thus  a  fixed  topography,  o  may  still  be  prescribed  arbitrarily 
and  k  and  £  determined  from 


k2  =  (a-1 )  0“'  /gH^  (a-2 ) 


(28) 


-a- 


and 

i.2  =  k2/ (oc-1)  =  a2/gHM(u-2)  (29) 

These  equations  are  analogous  to  eqs.  (8)  and  (9)  for  the  first  example, 
with  a  playing  the  role  which  played  there;  Che  -slope  at  the  coastline  in 
this  example,  ,  is  a  function  only  of  a. 


II :  Discussion  of  Solution 

The  topography  obtained  as  eq .  (19)  must  be  examined  to  determine  under 
what  conditions  it  would  be  physically  acceptable.  One  tool  to  this  end  is 
the  asymptotic  form  for  large  |£/HCD|,  e.g.  Abramowitz  and  Stegun  (1964): 
M(l,l+ct,z)  -(a/z)  (l-(l-a)/z  +  .  .);  z  <  0 


^  f  (l-ra)exp  (z)  (z)-a;  z  >  0 


(30) 


which  leads  to 


H(O/H0O  %  l  +  2  /  [  ( 1+v)  (£/Hod)  ]  +  .  .  ;  5/Hoo  <  0 

(31) 

^  -2F (a) exp (q/H  )(C/Hro)  +  .  .  ;  >  0 

Since  x  approaching  +<“  defines  the  appropriate  direction  for  the  ocean 
(a  decaying  edge  wave),  the  relationship  of  (  to  x  must  be  used  to  determine 
the  sense  of  F,. 

While  there  are  a  number  of  different  cases  to  be  discussed,  it  may  be 
easiest  to  first  consider  the  special  case  of  no  rotation,  f=0.  The  parameters 
a  and  b  will  be  taken  to  be  zero  throughout  for  convenience;  their  effect 
could  be  included  without  difficulty.  In  tiiis  case,  eq .  (21)  reduces  to 

ix  =  1  -  (f,/iU/(k2/i2  -  1)  =  1  -  (C/Hc) / (a-2 ) .  (32) 

For  a  >  2,  k2  must  be  greater  than  i2  and  is  positive;  the  sense  of  £  will 
then  be  oppostic  to  that  for  x.  For  1  u  <  2,  k2  must  be  less  than  £2  and  11*, 
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is  negative  (although  is  considered  as  an  asymptote  in  some  situations,  in 
reality  it  is  some  convenient  combination  of  symbols  which  measures  a 
horizontal  length  scale);  in  this  case,  the  sense  of  £  will  be  opposite  to  that 
for  x  but  Che  sense  of  (  £  / )  will  be  the  same  as  that  for  x.  Thus  for 
a  >  2,  eq .  (31)  yields  a  positive  depth,  H(£),  which  approaches  from  "below" 
(shallower  depths)  as  £  approaches  -<».  For  1  <  a  <  2,  eq .  (31)  yields  a 
positive  depth,  11(C),  which  increases  exponentially  as  £  approaches 
There  are  no  solutions  for  a  <  1  which,  have  real  f,  and  k.  For  a  =  2,  k2  equals 
£2  and  the  solution  for  11  must  be  redone,  leading  to  M(x)  =  (o2/g) (ix3/3  -  x2). 
Figure  1  sketches  the  nond imensional  depth,  H,  and  edge  wave  height,  E,  as  a 
function  of  £  for  a  =  4.0.  This  figure  is  the  same  for  all  frequencies; 
however,  when  a  dimensional  horizontal  length  scale,  £,  is  used,  the 
dimensional  depth,  H,  is  unchanged  but  the  edge  wave  height,  E,  depends  on 
frequency.  Figure  2  indicates  this  for  iUc.  =  200  meters,  a  =  4.0  with  E 
plotted  or.  a  relative  scale  for  two  periods,  43  sec.  and  i  sec . 

Finally,  Figure  3  shows  the  situation  when  =  -200  meters  and  a  =  1.5  for  a 
wave  period  of  43  sec.  This  topography  is  concave  downward,  i.e.  is  "dual" 
to  the  previous  topography.  Other  wave  periods  yield  different  curves  for 
E(£)  on  the  same  base  topography,  H(£). 

The  inclusion  of  rotation  introduces  some  new  effects  as  well  as 
modifications  of  the  non-rotating  solutions.  Again,  the  two  families 
(Ho,  >  0  and  Hcc  <  0)  are  considered  separately.  Before  proceeding,  it  is 
useful  to  provide  some  relationships  between  the  parameters,  e.g. : 
k'/P.2  =  v'-  +  (a-l)(l+v)2  =  a(l+v)2  -  1  -  2v 

a2/f2  =  ghOTja(l+v)2  -  1  -  2vj[.t(l+v)2  -  2  -  2v] / (a-1 ) (1+v) 2 
(k2/i2)(l-f2/,J2)  =  (,t-l)(l+v)2 
tx  =  (1+V)'1  - 


(C/lU)/[u(l+v)-2] 
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Consider  first,  the  conc.i'”  ..t,w.ird  family  corresponding,  Co  positive 
H^.  For  trequeucy  waves,  a  >  f ,  this  implies  k2  >  l2  which  can  take  on 

two  forms.  The  case  0  '  <  k  requires  v  >  0  and  a  >  2/(l+v)  and  yields  H(£) 

approaching  If,,  from  below  as  f  approaches  -<*  (i.e.  x  approaches  +**) . 

The  case  k  '  -1  <  0  requires  -1  •-  v  -  G  and  a  >  2/  (1+v) ;  a  >  1.  Again  H(0 

approaches  H  1  rom  below  as  (x)  approaciies  -(+)“. 

Next  consider  low  frequency  waves,  a  <  f.  This  implies  k2  <■  „2. 

For  0  <  k  i,  v  >  0  and  a  <  1,  (l+2v)/ (1+v)2  <  i  <  2/ (1+v)  to  have  real  k  and 
i.  Here  x,q  both  approacii  +«=  and  the  topography  is  negative  and  thus 
unacceptable.  Finally,  take  -l  <  k  <  0  which  requires  v  <  0  and  a  <  1, 

(l+2v)/ (1+v) z  <  a  <  2/(1 +v).  For  -1  <  v  <  0,  t,  approaches  +°°  as  x  approaches 

+00  and  a  negative  topography  is  again  obtained.  However  for  v  <  -1,  a  <  0  and 
4  approaciies  as  x  approaciies  +«.  Tiiis  provides  a  positive  topography,  but 
one  which  approaches  HOJ  from  above.  At  f  =  0,  however,  M(l,l+u,0)  =  1  and 
H(C=0)  =  I!a,i.  1  -  2/a  (1+v) ;  while  as  €,  approaciies  +■*•,  H(£)  is  negative.  Then  the 
topography  for  this  case  develops  a  trough  below  H  before  approaching  H 
from  above  (deeper  values). 

These  results  are  described  in  Figures  4,  5  and  6.  For  =  200  meters, 
a  =  2.5,  v  =  0.1,  Figure  4  shows  a  concave  upward  topography  similar  to  that 
of  Figure  2  as  does  Figure  5  for  =  200  meters,  a  =  2.5  and  v  =  -0.1. 

Figure  6  on  the  other  hand,  for  =  200  meters,  a  =  -0.77  and  v  =  -4.0 
indicates  a  new  kind  of  topography,  i.e.  one  with  a  trough.  These  figures  are 
shown  in  dimensional  variables  for  a  wave  period  of  48  sec,  but  could  be 
readily  scaled  to  other  periods  by  "stretching"  the  E  curves.  Figures  4  and  5 
contain  no  apparent  differences  from  tin;  results  shown  in  Figure  1;  in  all  of 
these,  the  edge  wave  decays  over  a  reasonable  portion  of  the  topography  and 


A 
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has  a  node  near  the  coastline.  Figure  6,  however,  is  cuite  different  in  that 
the  topography  is  not  monotonic,  but  contains  a  trough.  The  edge  wave  is 
concentrated  along  the  coastal  side  of  the  trough,  but  has  non-negligible 
values  on  the  other  side  as  well.  It  would  be  of  interest  to  make  a  paramet¬ 
ric  study  of  this  case  in  particular  since  analytical  solutions  for  non¬ 
monotonic  topographies  do  not  appear  to  currently  exist  in  the  literature. 

This  is  a  low  frequency  wave  traveling  "down"  the  coast  (k  <  0)  and  may  be 
related  to  the  trench  wave  of  Mysak,  LeBlond  and  Emery  (1979);  they  can  not 
be  the  same  since  tais  is  a  surface  gravity  wave  and  the  results  of  Mysak, 
et  al.  refer  to  non-divergent  (rigid  lid)  trapped  waves. 

N’ext ,  t  iie  i:-..,ily  of  convex  upward  topographies  corresponding  to 
negative  K^,  must  he  examined.  High  frequency  waves,  a  >  f,  require  k2  <  i.2. 

The  case  0  v  k  •'  i  requires  0  <  v  <  1  and  to  have  real  c  and  i  with 
a  >  1  requires  (l+2v)/ (1+v)2  <  a  <  2/(l+v).  As  x  approaches  +“,  (C/H*,) 
approaches  +*»  and  H ( q ) / is  negative  implying  H(£)  is  positive  and  increases 

exponentially  with  U/H^)  as  q  approaches  -<*>.  The  case  -SL  <  k  <  0  requires 

-1  <  v  <  0,  (l+2v)/ (1+v)2  <  a  <  2/(l+v)  and  a  >  1.  As  x  approaches  +®, 

(C/Hoo)  approaches  -K«  and  again  H(ti)  is  positive  and  increases  exponentially 
with  (C/H^)  as  £  approaches  -“>.  In  both  of  these  cases,  H(£=0)  is  negative. 

The  low  frequency  solutions,  a  <  f,  require  k2  >  l2.  For  0  <  SL  <  k, 
v  >  1  and  2/ (1+v)  <  a  <  1.  This  leads  to  (C/H^)  approaching  -»  as  x  approaches 
+°°  and  thus  to  an  unacceptable  negative  topography,  H(£).  Finally,  the  case 
k  <  -&  <  0  requires  v  <  -1  and  2/ (1+v)  <  a  <  1.  (£/Hm)  then  approaches  +»  as 

x  approaches  +°°  and  H(£)/HOT  approaches  sgn [+F (a) ]®.  Then  H(£)  will  be  positive 
(and  exponentially  increasing  with  -£)  if  T (a)  <  0  which  occurs  for  -1  <  a  <  0, 

-3  <  a  <  -2,  etc.  To  keep  1  +  a  >  0  (and  thus  M(l,l+a,z)  monotonic  as  z 
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approaches  +**),  only  the  case  -1  <  a  <  0  will  be  chosen.  Figures  7  and  8 
indicate  results,  in  nondimensional  form,  for  =  -200  meters,  a  =  1.1  and 
v  =  +0.1  and  -0.1  respectively.  These  results  indicate  an  edge  wave  which 
persists  over  a  significant  portion  of  the  topography  and  contains  a  node  near 
the  coastline.  Finally,  Figure  9  shows  results  for  =  -200  meters, 
a  =  -0.1  and  v  =  -1.1,  i.e.  the  last  case.  Here,  the  node  is  far  from  the 
coastline  and  the  edge  wave  is,  where  it  has  non-negligible  values,  negative 
(i.e.  out  of  phase  compared  to  the  other  figures)  over  most  of  the  topography. 


Ill:  Conclusion 

An  inverse  approach  to  long  edge  wave  trapping  has  resulted  in  a  new 
analytical  solution  as  well  as  an  alternative  approach  to  the  well-known 
Ball  (n=0)  solution.  These  new  solutions  appear  to  form  only  the  lowest 
(or,  in  view  of  the  single  node,  possibly  the  first)  modes  of  what  are  in 
reality  a  spectrum  of  edge  wave  modes;  this  may  be  seen  by  analogy  to  the 
Ball  solution,  and  a  reasonable  next  step  in  this  research  would  be  to 
extend  these  solutions  to  higher  modes. 

An  equally  reasonable  next  step  would  be  to  relate  these  mathematical 
solutions  to  physical  observations  to  determine  any  value  which  they  might 
have  beyond  their  novelty  as  an  "inverse"  view  of  the  world. 

Both  of  these  tasks  are  left  for  future  study  so  that  the  above 
results  may  at  this  time  at  least  be  presented  for  discussion. 


-  n  - 
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Figures 

1:  Nondimensional  depth  H  ( - )  and  wave  height  E  ( - )  as  functions 

of  £  for  v  =  0.0,  a  =  4.0  and  H  =  200  meters. 

oo 

2:  Dimensional  depth  H  ( - )  and  wave  height  E( - )  as  functions  of  £ 

for  v  =  0.0,  a  =  4.0,  =  200  meters  and  wave  periods  of  48  and  162 

seconds.  Note  that  H  is  independent  of  period  although  E  does  change 
(on  a  relative  scale) . 

3:  Dimensional  depth  H( - )  and  wave  height  E( - )  as  functions  of  £ 

for  v  =  0.0,  a  =  1.5  and  H  =  -200  meters  with  a  wave  period  of 
48  seconds.  This  topography  is  exponentially  increasing  with  distance 
from  shore,  but  could  be  terminated  at  a  flat  bottom  beyond  £  =  -1.5  km 
with  no  real  change.  Other  periods  would  yield  the  same  H  but  differing 
E. 


4:  Dimensional  depth  H  ( - )  and  wave  height  E  ( - )  for  the  rotational 

case  with  v  =  0.1,  a  =  2.5,  H  =  200  meters  and  a  period  of  43  seconds, 
k  =  0.44  x  10-4  cm-1,  Z  =  0.32  x  10-4  cm-1. 

5:  Dimensional  depth  H  ( - )  and  wave  height  E  ( - )  for  the  rotational 

case  with  v  «  -0.1,  a  =  2.5,  H  =  200  meters  and  a  period  of  48  seconds. 

'  '  OO 


Here  k  =  -0.68  x  10-4  cm-1  and  l  =  0.62  x  10~4 


6:  Dimensional  depth  H  ( - )  and  wave  height  E  ( - )  ior  v  =  -4.0, 

a  =  -0.77,  H  =  200  meters  and  a  period  of  48  seconds.  There  is  the 
"trough"  tonography :  note  that  E  is  still  significant  past  the  turning 
point  of  the  trough,  indicating  that  this  topography  is  "felt"  past  this 
point.  Note  that  k  here  is  negative,  i.e.  the  trapped  wave  progresses 
in  the  negative  (v)  coastline  direction  with  no  positive  k  counterpart 
as  was  the  case  for  Figures  4  and  5. 

7:  Nondimensional  depth  H  ( - )  and  wave  height  E  ( - )  for  v  =  0.1, 

a  =  1.1  and  =  -200  meters.  Again,  the  same  curves  hold  for  all  wave 
periods;  only  when  f,  is  dimensionalized  will  E  change  with  wave  period, 
k  here  is  positive. 

8:  Nondimensional  depth  H  ( - — )  and  wave  height  E  ( - )  for  v  =  -0.1, 

a  =  1.1  and  =  -200  meters.  This  figure  repeats  Figure  7  with  the 
sign  of  v  reversed.  The  results  are  quite  similar  although  here  k  will 
be  negative,  i.e.  a  counterpart  to  Figure  7. 

9:  Nondimensional  H  ( - )  and  E  ( - )  for  v  =  -1.1,  a  =  -0.1  and 

=  -200  meters.  Again,  k  is  negative  with  no  counterpart  for  positive 
k.  The  edge  wave  node  is  far  from  the  coastline  and  occurs  at  extremely 
small  (relative)  values  of  wave  height  giving  the  appearance  of  a 


monotonic  wave  height. 


200 


-0.2 


-4.0 


00= 200m 
v  =  -  4.0 
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